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Exact-diagonalization calculations for TV = 3 electrons in anisotropic quantum dots, covering a 
broad range of confinement anisotropies and strength of inter-electron repulsion, are presented for 
zero and low magnetic fields. The excitation spectra are analyzed as a function of the strength of the 
magnetic field and for increasing quantum-dot anisotropy. Analysis of the intrinsic structure of the 
many-body wave functions through spin-resolved two-point correlations reveals that the electrons 
tend to localize forming Wigner molecules. For certain ranges of dot parameters (mainly at strong 
anisotropy), the Wigner molecules acquire a linear geometry, and the associated wave functions 
with a spin projection S z = 1/2 are similar to the representative class of strongly entangled states 
referred to as VK-states. For other ranges of parameters (mainly at intermediate anisotropy), the 
Wigner molecules exhibit a more complex structure consisting of two mirror isosceles triangles. This 
latter structure can be viewed as an embryonic unit of a zig-zag Wigner crystal in quantum wires. 
The degree of entanglement in three-electron quantum dots can be quantified through the use of 
the von Neumann entropy. 
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I. INTRODUCTION 

Three-electron quantum dots are expected to attract 
a lot of attention in the near future due to several de- 
velopments, both experiment al an d theoretical. First, 
it was recently demonstratecP^El that detailed excita- 
tion spectra of two-electron quantum do ts (in addition 
to earlier ground-state measurements^^) can be mea- 
sured, and theoretically understood, as a function of the 
externally applied magnetic field. Thus, exploration of 
the excitation spectra of three-electron quantum dots ap- 
pears to be a next step to be taken. Second, three-qubit 
electron spin dev ices are expected to exhibit enhanced 
efficiencySEEEHo] f or quantum-computing and quantum- 
information purposes compared to single-qubit and two- 
qubit ones. 

In this paper, we carry out exact diagonalization 
(EXD) studies for a three-electron single quantum dot 
under low and moder ate magnetic fields. Unlike previ- 
ous EXD studies^J^ that focused mainly on the ground 
states (GSs) of circular quantum dotsp^l we investigate, 
in addition, the excitation spectra for three electrons in 
quantum dots with a wide range of anisotropies. More- 
over, consideration of anisotropic quantum dots allows us 
to investigate the structure of the many-body wave func- 
tions with respect to strong-correlations effects, such as 
electron localization and formation of Wigner molecules 
with a linear or zig-zag geometry. 

Most importantly, we investigate here the feasibility 
of producing mo del qua ntum entangled states (i.e., the 
socalled W stateipE^ES) ^ which are often employed in 
the mathematical treatment of quantum information and 
which have been experimentally realized with ultracold 
atoms in linear ion traps. 16 We note that a main factor 
motivating our investigations is the different nature of 
the entangling agent, namely, the electromagnetic field 



in the case of heavy ions versus the two-body Coulomb 
interaction in the case of electrons. 

We further mention other recent proposals in the con- 
text of solid state electronic devices for producing three- 
qubit entanglement. In particular, a scheme based on 
non-interacting electron-hole excitations in the Fermi sea 
was investigated in Ref. 17 Unlike our present study 
that focuses on the effect of the interparticle interaction, 
however, such interaction-free entanglement cannotP^ re- 
produce the symmetric VF-state [see Eq. (i"3| below]. A 
different proposal for realizing interaction-free entangle- 
ment uses pair-correlation functions to study tripartite 
entanglement shared among the spins of three fcrmions 
in a Fermi gas. 

The exact diagonalization method that we use for the 
solution of the Schrodinger equation corresponding to the 
Hamiltonian of three electrons interacting via a Coulomb 
potential in an anisotropic quantum dot, in conjunction 
with an analysis employing spin-resolved two-point cor- 
relation functions, allows us to gain deep insights into the 
nature of electronic states and three-qubit entanglement 
in real solid state devices. Additionally, the EXD method 
provides benchmark results, which could be used for as- 
sessment of the adequacy and relative accuracy of certain 
approximation schemes, including the model Heisenberg 
Hamiltonian for three localized electrons arranged in a 
ring geometry that was most recently used in an inves- 
tigation of the entangled ground states in a thrce-spin- 
qubit system.^ 



II. OUTLINE OF THE EXACT 
DIAGONALIZATION MANY-BODY METHOD 

We consider three electrons under zero or low magnetic 
field (B) in a single quantum dot. The corresponding 
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many-body Hamiltonian is written as 
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where k is the dielectric constant of the semiconductor 
material (12.5 for GaAs). The single-particle Hamilto- 
nian is given by 



H = T+V(x,y)+g* f i B Ba, 



(2) 



where the last term is the Zeeman interaction, with g* 
being the effective Lande factor, \xb the Bohr magneton, 
B the perpendicular magnetic field, and a — ±1/2 the 
spin projection of an individual electron. The kinetic 
contribution in Eq. Q is given by 



T 



(e/c)A(r)]- 
2m* 



(3) 



with to* being the effective mass (0.067m e for GaAs) 
and the vector potential A(r) = 0.5(— Byt + Bxj) being 
taken according to the symmetric gauge. The external 
confining potential is denoted as V(x, y), where r = xi + 

The external potential is modeled by an anisotropic 2D 
oscillator 
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V(x,y) = ~™*(u x x +lu y ), 



(4) 



which reduces to a circular parabolic QD for u) x = u) v = 
u)q. The ratio rj = uj x /lo v characterizes the degree of 
anisotropy of the quantum dot, and it will be referred 
to thereafter as the anisotropy parameter. Results will 
be presented for three cases: rj = 1 (circular), rj = 0.724 
(slightly anisotropic), and rj =1/2 (strongly anisotropic). 

We find the eigenstates of the many-body Hamiltonian 
([lj using an exact diagonalization method. Accordingly, 
we expand the many-body wave function as a linear su- 
perposition, 



EXD 



(ri,r 2 ,r 3 ) 



l<i<j<k<2K 

(5) 

where \ip(l;i)ip(2;j)tp(3;k)) denotes a Slater determi- 
nant made out of the three spin-orbitals t/j(l;i), tp(2;j), 
and tp(3;k). For the spin orbitals, we use the notation 
V>(1; t) = v>i(l T) if 1 < i < K and V(l; t) = Vi-jf (1 i) if 
K+l <i<2K [and similarly for ip(2;j) and ^(3; k)]. K 
is the maximum number of space orbitals <Pi(v) that are 
considered, with ipi{l |) = ipi(ri)a and ipi(l I) = ipi(r{)l3 
where a and (3 denote up and down spins, respectively. 
The space orbitals <pi(r) are taken to coincide with the 
real cigcnfunctions of a 2D anisotropic oscillator, that 
is, the index i = (m,n) and <^i(r) = X m (x)Y„(y), with 
X m (Y n ) being the cigcnfunctions of the corresponding 
one-dimensional oscillators in the x{y) direction with fre- 
quency uj x {uj y ). The parity operator V yields VX m (x) = 
( — l) m X m (x), and similarly for Y n (y). 
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FIG. 1: (Color online) Ground-state and excitation en- 
ergy spectra [referenced to 3ft \/ti;Q + wf/4, with ojq = 
■sj {uj% + ujy)/2] as a function of the magnetic field for N = 3 
non-interacting electrons in a circular quantum dot (r) = 1). 
Parameters: external confinement huj x = hui y = 5 meV; di- 
electric constant k — oo; effective mass m* = 0.067m e , ef- 
fective Lande coefficient g* — 0. The labels (S; L) denote 
the quantum numbers for the total spin and the total angular 
momentum. Different Landau bands are denoted by the dif- 
ferent M values. The S z indices are not indicated, since the 
multiplets (S = 1/2, Sz) and (5 = 3/2,52) are degenerate in 
energy when g* = 0. 



The total energies -Eexd and the coefficients Aijk's 
are obtained through a direct numerical diagonalization 
of the matrix eigenvalue equation corresponding to the 
Hamiltonian in Eq. ([T]). For the solution of this large 
scale, but sparse, matrix eigenvalue problem, we have 
used the ARPACK computer codePl 

The EXD wave function ^ preserves by construction 
the third projection S z of the total spin, since only Slater 
determinants with a given S z value are used in the ex- 
pansion. The exact diagonalization automatically pro- 
duces eugenfunctions of the square, S 2 , of the total spin 
S = J2i=i^i- The corresponding eigenvalues S(S + 1) 
are calculated with the help of the expression 



S 2 |SD) = 



(N a - iV^) 2 /4 + JV/2 + £- 



i<j 



|SD), (6) 



where |SD) denotes a Slater determinant and the oper- 
ator Wij interchanges the spins of electrons i and j pro- 
vided that their spins are different; N a and Np denote 
the number of spin-up and spin-down electrons, respec- 
tively, while ./V denotes the total number of electrons. 

Since the spin orbitals ijj's are orthogonal, the Coulomb 
matrix elements between two Slater determinants are cal- 
culated using the Slater rulesj^ and the necessary two- 
body matrix elements between space orbitals 



J J dridr 2 ^*(ri)^*(r 2 )^-^- 



r 2 



r (pfc(ri)^(r 2 ) (7) 
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are calculated numerically. We have found that this 
method produces numerically stable results in compar- 
ison with algebraic expressions.^ 



III. ENERGY SPECTRA 

In this section, we study the ground-state and excita- 
tion spectra as a function of an increasing magnetic field 
B with an emphasis on the role of correlation effects and 
the influence of the anisotropy. 

To better understand the importance of correlations, 
we first display in Fig.[l]the spectra in the absence of the 
Coulomb interaction (non-interacting electrons) and for 
the case of a circular quantum dot. These energy spectra 
can be determined simply as X)i=i e ? F (^)> where e^ F (B) 
are the Darwin-Fock energies for a single electron ) 21 * 22 * 23 ! 
The main trend is the formation of three-particle Landau 
bands (each with an infinite number of states) that tend 
for B — > oo to the asymptotic energy levels (A4+3/2)hui c , 
Ai = 0,1,2,.... Note that, for large magnetic fields 
(B — ► oo), the reference energy 3hy/u)Q + w 2 /4, with 

luq = + ujy)/2], tends to 3Hcu c /2. In this limit, 

the states (S,L), belonging to the same Landau band 
M, become degenerate in energy, converging to the cor- 
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FIG. 2: (Color online) Ground-state and excitation en- 
ergy spectra [referenced to 3ft ^/(Jq + tjf/4, with uiq = 
sj (co'i + lj%)/2] as a function of the magnetic field for N = 3 
non-interacting electrons in a anisotropic quantum dot with 
anisotropy parameter (77 — 0.724). Parameters: external con- 
finement fiui x — 4.23 meV; hu y = 5.84 meV; dielectric con- 
stant k = 00; effective mass m* = 0.070m e , effective Lande 
coefficient g* = 0. The labels (S; L) denote the quantum 
numbers for the total spin and the total angular momentum 
in the corresponding circular quantum dot. Although the to- 
tal angular momentum is not a good quantum number for an 
anisotropic quantum dot, we retain the labels L here in order 
to facilitate comparison with the circular case in Fig. [I] The 
S z indices are not indicated, since the multiplets (S = 1/2, S z ) 
and (S = 3/2, S z ) are degenerate in energy when g* — 0. 
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FIG. 3: (Color online) Ground-state and excitation en- 
ergy spectra [referenced to 3ft^/wQ + uf/4, with uiq = 
*/ (uj\ + ujy)/2] as a function of the magnetic field for TV = 3 
interacting electrons in a circular quantum dot (j) = 1). Pa- 
rameters: external confinement hui x = hw y = 5 meV; dielec- 
tric constant k — 12.5; effective mass m* = 0.067m e ; effective 
Lande coefficient g* = 0. The labels (S; L) denote the quan- 
tum numbers for the total spin and the total angular momen- 
tum. The S z indices are not indicated, since the multiplets 
(S — 1/2, S z ) and (S = 3/2, S z ) are degenerate in energy 
when g* = 0. 



responding familiar Landau level (with index M.). Apart 
from an overall constant, the picture in Fig. [l] is the 
same as that found in the phcnomcnological "constant- 
interaction" modelP^ An important property is the ab- 
sence of crossings between individual levels within each 
Landau band. A consequence of this is that the ground 
state at any B has the same quantum numbers as the 
one at B = 0, i.e., it has total spin S — 1/2 and total 
angular momentum L = 1. 

The absence of crossings within each Landau band 
is a characteristic property of non-interacting electrons, 
and it is independent of the anisotropy of the external 
confinement. This point is illustrated in Fig. [2] where 
the non-interacting three-electron spectra are plotted for 
the case of a quantum-dot with moderate anisotropy 
(?/ = 0.724). [for the single-electron energies ef F (B) in 
an elliptic quantum dot, see Ref . EMI and Ref . 125] . An in- 
spection of Fig. [2] shows that the anisotropy has an effect 
mainly for small values of the magnetic field (by lifting 
the degeneracies at B — 0). On the formation of the Lan- 
dau bands at higher B, the anisotropy has practically no 
effect, and, in particular, it cannot induce level crossings 
within each Landau band) . 

Another property of the non-interacting spectra is the 
existence of several degenerate levels (not shown in Figs, 
[l] and [2]) associated with the excited states. We have 
searched for such degeneracies by inducing a small lifting 
of them through the artificial use of a very weak Coulomb 
repulsion specified by k = 200. For example, in the cir- 
cular case [see Fig. [I], we found that the state (1/2; 2) is 
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doubly degenerate, while the state (3/2; 3) is degenerate 
with two other (1/2; 3) states. These additional states 
move higher in energy as the strength of the Coulomb 
interaction increases. We further found that the lifting 
of degeneracies is sufficiently strong for larger Coulomb 
repulsions with k < 12.5 that all the curves in Figs.[3j[4] 
and [5] below are simple (i.e., the additional states present 
in the non-interacting case have been pushed much higher 
and fall outside the energy window shown). 

Turning on the interaction introduces correlation ef- 
fects that lead to important modifications of the non- 
interacting spectra shown in Figs, [l] and [2] Fig.[3]displays 
the corresponding spectra for the same circular quantum 
dot as in Fig. [T] but in the presence of a Coulomb re- 
pulsion with k — 12.5 (GaAs). Of course, a first effect 
is the increase in the total energy, but the main differ- 
ence from the non-interacting case in Fig. [l] is the pres- 
ence of crossings between levels within the same Lan- 
dau band. As a result, within the plotted range of mag- 
netic fields, the ground-state total-spin quantum number 
remains 5 = 1/2 at the first ground-state crossing (at 
point B), and then it changes to 5 — 3/2 (at the second 
ground-state crossing at point C) . At the same time, the 
total angular momentum changes from L — 1, to L = 2, 
and then to L = 3, respectively. As long as the effective 
Lande coefficient g* = 0, which is the case for the results 
presented in this section, this threefold alternation in the 
spin and angular momentum quantum numbers repeats 
itself ad-infinitum. We note that experimental observa- 
tion of this threefold alternation may be forthcoming, 
since quantum dots with a vanishing Lande coefficient 
have been recently fabricated^ and were used already to 
measure two-electron excitation spectra. 

The crossings of the curves associated with the three 
different pairs of quantum numbers (5 = 1/2; L = 1), 
(1/2; 2), and (3/2; 3) form a small triangle (labeled as 
ABC), which is located about B ~ 3.4 T. Anticipating 
the results for non-circular dots below, we note that this 
triangle tends to collapse to a single point with increasing 
anisotropy. 

Another prominent difference between the spectra of 
non-interacting (Fig. [IJ and interacting (Fig. [3]) elec- 
trons pertains to the degeneracies at B — between 
the 5 — 3/2 and S — 1/2 states that are lifted in 
the interacting-electrons case; compare in particular the 
curves with quantum numbers (1/2; 2) and (1/2; 0) with 
the (3/2; 0) one. In contrast, the original degeneracies at 
B = of the S — 1/2 states are unaffected by the in- 
terelectron interaction; compare the curves (1/2,1) and 
(1/2, -1), as well as the ones labeled (1/2, 2) and (1/2, 0). 
However, these S — 1/2 degenaracies at B — are lifted 
as a result of an increasing anisotropy of the quantum 
dot, as seen in Fig. [3] 

Next, we explore the effect of increasing the anisotropy 
of the quantum dot. In particular, keeping the same 
strength for the Coulomb interaction (n = 12.5), we 
present two representative anisotropy cases, i.e., r\ = 
0.724 (intermediate anisotropy, see Fig. H, and 77 = 1/2 



(strong anisotropy closer to a quasilinear case, see Fig. 

Inspection of the results for the case of intermediate 
anisotropy (Fig. [4j , reveals that compared to Fig. [3] the 
spectra are distorted, but they maintain the overall topol- 
ogy of the circular dot. As a result, we have been able 
to use the same pairs of labels in naming the different 
curves, even though the second label does not have the 
meaning of an angular momentum (the total angular mo- 
mentum is not conserved for rj =/= 1). There are two main 
differences from the circular case: i) the degeneracies at 
B = between the S =1/2 states are lifted, and ii) there 
is a marked rounding of all the 5=1/2 curves in the be- 
ginning, so that they do not intersect the vertical energy 
axis at sharp angles as is the case with Fig. [3j This ini- 
tial rounding and bending of the energy curves due to 
the anisotropy has been experimentally observed^SI m 
two-electron quantum dots. 

In the case of strong anisotropy (Fig. [5|, the spectra 
have evolved to such an extent that only little relation 
to the circular case can be traced, and as a result we 
use a single label signifying the total spin to distinguish 
them. An important feature that emerges is that the 
three curves with lowest energies (two 5=1/2 and one 
5 = 3/2 curve) form a band that is well separated from 
the other excited states. The existence of such an iso- 
lated lowest-energy band is important for validating sim- 
ple two-qubit and three-qubit models introduced in quan- 
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FIG. 4: (Color online) Ground-state and excitation en- 
ergy spectra [referenced to Shy^uig + u/iji, with ujo = 
v/ (w| + uiy)/2] as a function of the magnetic field for TV = 3 
interacting electrons in an elliptic quantum dot with interme- 
diate anisotropy (anisotropy parameter 77 = 0.724). Parame- 
ters: external confinement hu) x = 4.23 meV, huj y = 5.84 meV; 
dielectric constant k = 12.5; effective mass m* = 0.070m e ; ef- 
fective Lande coefficient g* = 0. The labels (S; L) denote the 
quantum numbers for the total spin and the total angular mo- 
mentum in the corresponding circular quantum dot. The S z 
indices are not indicated, since the multiplets (S = 1/2, S z ) 
and (S = 3/2, S z ) are degenerate in energy when g* = 0. Note 
the shrinking of the ABC triangle compared to the f] = 1 case 
shown in Fig. [3] 
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FIG. 5: (Color online) Ground-state and excitation en- 
ergy spectra [referenced to 3ft -^/cjq + wf/4, with ljq = 
■sj (uj% + ujy)/2] as a function of the magnetic field for N = 3 
electrons in an elliptic quantum dot with strong anisotropy 
(anisotropy parameter r\ = 1/2). Parameters: external con- 
finement hu) x = 3.137 meV, hoj y = 6.274 meV; dielectric con- 
stant re = 12.5; effective mass m* — 0.067m e ; effective Lande 
coefficient g* — 0. The single labels denote the quantum num- 
bers for the total spin. The S z indices are not indicated, since 
the multiplets (S = 1/2, S z ) and (S = 3/2, S z ) are degenerate 
in energy when g* = 0. Note the collapse of the triangle ABC 
compared to the cases with 7] = 1 (Fig. [3} and r) = 0.724 (Fig. 
E|, and the appearance of a triple-point crossing. 



turn computation and quantum information theoryP^ll 
Another remarkable feature of the strong-anisotropy 
case is the appearance of a non-trivial triple-point cross- 
ing lying on the ground-state curve (see arrow in Fig. 
[5|, which is created from the collapse of the ABC tri- 
angle between the two S = 1/2 and the one S = 3/2 
lowest-in-energy curves (compare Figs. |3]and|4|. This 
low-energy non-trivial triple point [forming within the 
lowest Landau band (Ai = 0)] is due to the effect of the 
Coulomb interaction, and it is to be contrasted to other 
trivial triple-point crossings at much higher energies aris- 
ing from the intersection of the lowest Landau band with 
the M. — 1 and M. = 2 higher Landau bands, and which 
are present even in the non-interacting limit [see, e.g., the 
triple crossing at (2.0 T; 14.2 meV) in Fig. [I]. It would 
be of interest to analyze whether the recently observed^ 
triple-point crossings in deformed quantum dots are non- 
trivial or trivial in the sense described above. 

Before leaving this Section, we note that the spin mul- 
tiplets (S — 1/2, S z ) and (S = 3/2, S z ) are degenerate in 
energy when g* = 0, which was the case for the energy 
spectra presented in Figs. 1-5. At a given magnetic field, 
this degeneracy is naturally lifted when g* ^ 0; how- 
ever, the final total energies can be easily calculated by 
adding the Zeeman term g* fisBS z to the spectral curves 
displayed in these figures. Furthermore, for a given pair 
(S, S z ), the Zeeman term does not influence the intrinsic 
structure of the many-body EXD wave function [i.e., the 
expansion over constituent Slater determinants, see Eq. 




FIG. 6: (Color online) Exact-diagonalization electron densi- 
ties for the ground states of N = 3 electrons in an anisotropic 
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dot with parameters tuo x = 3.137 meV, 
(r) = 1/2), effective mass m* = 0.067ra e , dielectric constant 
re = 12.5 (GaAs). (a): The case of zero magnetic field, B — 0. 
(b) The case with a magnetic field B — 6 T. Lengths in nm. 
The electron densities are in arbitrary units, but with the 
same scale in both panels. 



([5J], and thus taking g* — does not effect the results for 
electron densities, conditional probability distributions, 
and von Neumann entropies presented below. 



IV. MANY-BODY WAVE FUNCTIONS FOR 
STRONG ANISOTROPY (77 = 1/2) 

A. S = 1/2 ground states: Evolution of electron 
densities as a function of the inter-electron repulsion 

When the confining potential lacks circular symme- 
try, charge localization is reflected directly in the single- 
particle electron densities. Indeed, electron localization 
is visible in Figs. [6] and [7] which display the electron den- 
sities for N = 3 electrons in an anisotropic quantum dot 
with rj = 1/2. Fig. [6] illustrates the evolution of electron 
localization with increasing magnetic field in the case of 
a weaker Coulomb repulsion (k — 12.5). One sees that 
already at B = 0, the electron density is shaped linearly 
for all practical purposes. However, the three peaks of 
the localized electrons are rather weak, which contrasts 
with the case of B — 6 T [Fig. |6jb)], where the three 
electron peaks are sharply defined. 

Fig. [7] [in conjunction with Fig. [6]ja)] illustrates the 
strengthening of electron localization as a function of 
increasing Coulomb repulsion, i.e., decreasing dielectric 
constant k, from a value of 12.5 [Fig.^a)] to k = 3 [Fig. 
1 'a)] and then to k = 1 [Fig. [7^b)]. In this last case [Fig. 
7'b)], the three electrons arc almost fully localized, with 
orbitals that exhibit practically zero mutual overlap. 

Since we keep the average frequency, ojq = 
\J (w| + lu%) /2, approximately constant (i.e., Kuq ps 5.0 
meV) for all anisotropy cases studied in this paper, de- 
creasing the dielectric constant is equivalent to increas- 
ing the Wigner parameter^! R w . At zero magnetic field, 
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FIG. 7: (Color online) Exact-diagonalization electron densi- 
ties at zero magnetic field (B = 0) for the ground state of 
N = 3 electrons in an anisotropic dot with parameters hui x = 
3.137 meV, hu y = 6.274 meV (jj = 1/2), m* = 0.067m,;. (a): 
dielectric constant k = 3.0. (b): dielectric constant k — 1.0. 
Lengths in nm. The electron densities are in arbitrary units, 
but with the same scale as in Fig. [6] for both panels. 



Rw is widely used as a universal parameter to indicate 
the strength of correlations, since it provides the relative 
strength of the Coulomb repulsion with respect to the 
quantum kinetic energy, i.e., 



Rw = 



7(^ ) 

Hluq 



(8) 



with the characteristic length Iq = \Jhj (m*uo). For the 
numerical values of Rw associated with the cases studied 



here, see Fig. 15 



B. 5 = 1/2 ground state: Spin resolved intrinsic 
structure for strong repulsion (re = 1) 

In the previous section, we saw that already the elec- 
tron densities provide partial information about the for- 
mation of a linear Wigncr molecule within an elliptic 
quantum dot. Indeed, from the charge distributions in 
Figs. [6] and [7] one can infer that the electrons are localized 
in three separate positions Ri, R 2 , and R 3 . If the elec- 
trons were spinless, this situation could be approximately 
reproduced by a single Slater determinant denoted as 
| O OO)- However, to probe the spin distribution of 
the electrons, the exact-diagonalization charge densities 
do not suffice; one needs to consider spin-resolved two- 
point correlation functions, defined as 

P aao (r,r ) = 

(^exd| £ s{r _ ri)s(rQ _ rj)Saai S aoryj |* EXD ), (9) 

with the EXD many-body wave function given by equa- 
tion 

Using a normalization constant 



Af(a,a ,r ) = J P ffo - (r,r )<ir, 



we further define a related conditional probability distri- 
bution (CPD) as 

V aao (r,r ) = P CTCTo (r,r )/7V(cr,(T ,ro), (11) 

having the property J V aao (r, r )dr = 1. The spin- 
resolved CPD gives the spatial probability distribution 
of finding the remaining electrons with spin projection 
a under the condition that the third electron is located 
(fixed) at ro with spin projection do; cr and <tq can be 
either up (T) or down (|). 

Before examining such CPDs, evaluated for numeri- 
cally determined EXD wave functions, it is instructive 
to consider on a qualitative level the spin structure of 
the wave functions that can be formed from three local- 
ized spin-orbitals only. In particular, we focus on the 
case with a total spin projection S z = 1/2, when the 
most general three-orbital wave function is given by the 
superposition of three Slater determinants, i.e., by the 
expression 



$(£, = !) = o| UT) + &I TU) + c| ITT), 



(12) 



with the normalization a 2 + b 2 + c 2 = 1 . Unlike the circles 
used earlier to indicate spinless electrons, the arrows in 
Eq. ( 12 1 indicate the spin projections of the individual 



spin-orbitals. 



The general states ( 12 1 are a superposition of three 



Slater determinants and have attracted a lot of atten- 
tion in the mathematical theory of entanglement. In- 
deed, they represent a prototypical class of three-qubit 
entangled states known as W-states.^ For general coef- 
ficients a, b, and c, the states ( 12 ) are not eigenfunc- 
tions of the square of total spin S 2 (while the exact- 
diagonalization wave functions in Eq. ([5] are always good 
eigenfunctions of S 2 ). However, the special values of 
these coefficients that lead to g ood total-spin quantum 
numbers can be determined ) 7 * 30 ! In particular, using the 
notation $(S, S z ; i) (where the index i is employed in case 
of a degeneracy) , one has 

V3*(f,i) = |UT> + im> + UTT) (13) 

(i.e., a = b = c = l/y/3), 

V6$(|,|;l) = 2|UT>-|m>-UTT) (14) 
(i.e., a = 2/\/6, b = c = -1/V&), 

\/2<&(A,|;2) = |m>-|m) (15) 

(i.e., a = 0, b = l/y/2, c = -l/y/2). 

For completeness, we list the case for three fully spin- 
polarized localized electrons (which of course is not a 
W-state). 



*(§,§; 



TTT 



(16) 



(10) 



The wave functions with projections S z = —1/2 and S z = 
—3/2 are similar to the above, but with inverted single- 
particle spins. 
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FIG. 8: (Color online) Spin-resolved conditional probability 
distributions for the (1/2,1/2) ground state of N = 3 electrons 
in an anisotropic dot at zero magnetic field (B = 0) with pa- 
rameters huj x — 3.137 meV, fkj y = 6.274 meV (r; = 1/2), 
m* = 0.067m e and k = 1 [for the corresponding electron den- 
sity, see Fig. [7jb)] . The heavy arrow (green online) indicates 
the location of the fixed electron at ro [see Eq. |TTJ], with the 
indicated spin projection ao, i.e., up (|) or down (J.), (a) fj. 
CPD with the fixed spin-down electron located at the center, 
(b) 11 CPD with the fixed spin-down electron located on the 
right, (c) 11 CPD with the fixed spin-up electron located on 
the right, (d) J. j CPD with the fixed spin-up electron located 
on the right. The spin of the fixed electron is denoted by a 
thick arrow (green online). Lengths in nanometers. The ver- 
tical axes are in arbitrary units, but the scale is the same for 
all four panels. 



Before proceeding further, we note that the term Te- 
state is some times reserved for the symmetric form 
<&(§, |) in Eq. (13). This symmetric W-state has been 
experimentally realized in linear ion traps As we show 
below, quantum dots offer the means for generating in 
addition the less symmetric forms given by Eqs. ( 14 1 and 
(151. 



Such nonsymmctric three-qubit states are some- 
time denoted^ as W' -states (with a prime). In this pa- 
per, we do not make use of this distinction, and we refer 
to both symmetric and nonsymmetric forms simply as 
fF-states. 

In Fig. [8] we present several spin-resolved CPDs asso- 
ciated with the EXD ground state at B = and strong 
anisotropy r) = 1/2, which is a \I/ EXD (l/2, 1/2) state [see 




FIG. 9: (Color online) Spin-resolved conditional probability 
distributions for the (1/2,1/2) first excited state of N = 3 
electrons in an anisotropic dot at zero magnetic field (B = 0) 
with parameters huj x = 3.137 meV, hu y = 6.274 meV (-q = 
1/2), m* = 0.067m e and k = 1. (a) electron density (ED), 
(b) 11 CPD with the fixed spin-down electron located on the 
right at (61,0). (c) 11 CPD with the fixed spin-up electron 
located on the right at (61,0). (d) H CPD with the fixed 
spin-up electron located at the right at (61,0). The spin of 
the fixed electron is denoted by a thick arrow (blue online). 
Lengths in nanometers. The vertical axes are in arbitrary 
units, but the scale is the same for all four panels. 



Fig.|]. Although the EXD expansion in Eq. (5]) consists 
of a large number of Slater determinants built from delo- 
calizcd harmonic-oscillator orbitals, the CPD patterns in 
Fig. [8] reveal an intrinsic struct ure similar to that of the 
wave function <!>( | , | ; 1) in Eq. ( 14 1, which is made out of 
only three localized spin-orbitals. In particular, when one 
requires that the fixed electron has a down spin and is lo- 
cated at the center of the quantum dot, the spin-up elec- 
trons are located on the left and right with equal weights 
[Fig. |8ja)]. Keeping the down spin-direction, but moving 
the fixed electron to the right, reveals that the spin-up 
electrons are located on the left and the center with equal 
weights [Fig. pjb)]. Considering a spin- up direction for 
the fixed electron and placing it on the right reveals that 
the remaining spin-up electron is distributed on the left 
and the center of the quantum dot with unequal weights; 
approximately 4 (left) to 1 (center) following the square 
of the coefficients in front of the determinants | tit ) 
(a = 2/V6) and | J. Til (c = 1/V6) in the wave function 
$(2, 5)1) [ see Eq. ( 14 1] . Similarly, considering a spin- 
up direction for the fixed electron and placing it on the 
right reveals that the spin-down electron is distributed 
on the left and the center of the quantum dot with un- 




FIG. 10: (Color online) Spin-resolved conditional probability 
distributions for the (3/2,1/2) second excited state of N = 3 
electrons in an anisotropic dot at zero magnetic field (B — 0) 
with parameters hui x = 3.137 meV, hio y = 6.274 meV (77 = 
1/2), m* = 0.067m e and k = 1. (a) || CPD with the fixed 
spin-up electron located on the right at (70,0). (b) || CPD 
with the fixed spin-up electron located on the right at (70,0). 
(c) | j CPD with the fixed spin-down electron located on the 
right at (70,0). (d) || CPD with the fixed spin-down electron 
located at the center. The spin of the fixed electron is denoted 
by a thick arrow (green online). Lengths in nanometers. The 
vertical axes are in arbitrary units, but the scale is the same 
for all four panels. 



equal weights - approximately 1 (left) to 4 (center), in 
agreement with the weights of the Slater determinants in 
Eq. pi). 



C. 5=1/2 first excited state: Spin resolved 
intrinsic structure for strong repulsion (k = 1) 



In section |TV B| we investigated the intrinsic structure 
of the ground-state three-electron wave functions with 
total spin S =1/2 and for the case of a strong anisotropy 
i] = 1/2. In this section, we analyze a case of the first- 
excited EXD wave function with total spin S = 1/2 and 
for the same strong anisotropy r\ = 1/2, again at B = T 
and for strong interelectron repulsion k = 1. We denote 
this state as * EXD (l/2, 1/2; 2). 

In Fig. [9^ a), we display the electron density (ED) for 
this second S = 1/2 state, while in Figs. [9^b,c,d), we 
display spin-resolved CPDs for the same state. From 
the charge density, we conclude that the three electrons 
form a sharply defined linear Wigner molecule. The spin- 
resolved CPD with a spin-down fixed electron placed on 




FIG. 11: (Color online) Spin-resolved conditional probability 
distributions for the (3/2,1/2) ground state of N = 3 electrons 
in an anisotropic dot at B = 5 T with parameters huj x = 4.23 
meV, hujy = 5.84 meV (77 = 0.724), m* = 0.070m e and n = 
12.5. (a) || CPD with the fixed spin-up electron located 
on the right at (30,0). (b) || CPD with the fixed spin-up 
electron located on the right at (30,0). (c) || CPD with the 
fixed spin-down electron located on the right at (30,0). (d) || 
CPD with the fixed spin-down electron located at the center. 
The spin of the fixed electron is denoted by a thick arrow 
(blue online). Lengths in nanometers. The vertical axes are 
in arbitrary units, but the scale is the same for all four panels. 



the right [see Fig. [9[b)] is similar to that in Fig. [8j^ 
However, the two spin-resolved CPDs with a spin-up 
fixed electron placed on the right [see Figs.[9jc) and[9^d)] 
are quite different from the corresponding CPDs in Figs. 
|8jc) and [8jd). In fact, in both cases, only one single 
hump appears to the left of the fixed electron, located at 
the center for the remaining spin-up electrons [Fig. KHc)], 
or on the left for the remaining spin-down electrons [Fig. 
Id)]. 

This indicates that the intrinsic structure of the 
\I> EXD (l/2, 1/2; 2) wave function is close to that of 



$(§, §;2) in Eq. (15 1, with a = and b = -t 



D. S — 3/2 second excited state: Spin resolved 
intrinsic structure for strong repulsion (k = 1) 



In sections |l V B | and |I V C | we investigated the intrinsic 
structure of the many-body three-electron wave functions 
with total spin S = 1/2 and for the case of a strong 
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B=5T K=12.5 



B=0 K =1 




FIG. 12: (Color online) Spin-resolved conditional probability 
distributions for the (3/2,1/2) ground state of N — 3 electrons 
in an anisotropic dot at B = 5 T with parameters hui x = 4.23 
meV, huy = 5.84 meV (r] = 0.724), m* = 0.070m e and k = 
12.5. (a) || CPD with the fixed spin-down electron located 
on the y-axis at (0,20) (solid dot), (b) || CPD with the fixed 
spin-down electron located on the y-axis at (0,-20) (solid dot). 
The spin of the fixed electron is denoted by a thick arrow 
(blue online). Lengths in nanometers. The vertical axes are 
in arbitrary units, but the scale is the same for all panels in 
this figure and in Fig. |11| 



anisotropy r\ = 1/2. In this section, we analyze a case 
of an EXD wave function with total spin S — 3/2 and 
for the same strong anisotropy r\ — 1/2, again at B = 
T. In particular, we analyze the intrinsic structure of a 
^l/ EXD (3/2, 1/2) wave function that is the second excited 
state for these parameters. 

In Fig. [lOj we display spin-resolved CPDs for this 
5 = 3/2 excited state. A remarkable feature is that 
for a fixed electron placed on the right all three CPDS, 
TT [Fig. ©a)], It [Fig.^b)], and || [Fig. [Tofc)] coin- 
cide. This indicates that the intrinsic structure of the 
* EXD (3/2, 1/2) wave function is close to that of $(§, ~) 
in Eq. (U3|, with all three coefficients equal to each other, 
i.e., a = b = c. 



Taking into account the j| CPD with the fixed electron 
at the center of the quantum dot, it is clear that the 
geometric arrangement of the three localized electrons is 
linear. Arrangements that are more complicated than 
the linear geometry can emerge, however, for a range of 
different parameters, as is discussed in section [V A| below. 



V. MANY-BODY WAVE FUNCTIONS FOR 
INTERMEDIATE ANISOTROPY (rj = 0.724) 

A. Moderate repulsion (re = 12.5) 

In this section, we analyze a case of an EXD wave func- 
tion with total spin 5 = 3/2 and for the intermediate 
anisotropy rj — 0.724. In particular, we analyze the in- 
trinsic structure of a \l/ EXD (3/2, 1/2) wave function that 
is the ground state at a magnetic field B = 5T (see Fig. 



FIG. 13: (Color online) Exact-diagonalization electron densi- 
ties for the ground state of TV = 3 electrons in an anisotropic 
quantum dot with parameters hu x = 4.23 meV, hui y = 5.84 
meV (rj — 0.724, intermediate anisotropy) and m* = 0.070m e . 
(a) the (3/2,1/2) ground state at B = 5 T and re = 12.5. (b) 
the (1/2,1/2) ground state at 73 = and re = 1 (strong in- 
terelectron repulsion). Lengths in nm. The electron densities 
are in arbitrary units, with a different scale in each panel. 



In Fig. 11 we display spin-resolved CPDs for this 
ground state. A remarkable feature is that for a 
fixed electron placed on the right all three CPDS, |T 
[Fig. gja)], IT [Fig. [HJb)], and || [Fig. [TTJc)] coin- 
cide. This indicates that the intrinsic structure of the 
* EXD (3/2, 1/2) wave function is close to that of $(§, §) 
in Eq. ( fl3| ), with all three coefficients equal to each other, 
a = b = c. 

However, these CPDs, as well as the || CPD with 
the fixed spin-down electron at the center [Fig. [ITjd)], 
are broader along the y-direction compared to the CPDs 
associated with the linear molecular arrangement in Fig. 
1 1 0| This suggests that, for an intermediate anisotropy 
(77 = 0.724), the intrinsic structure of f EXD (3/2, 1/2) is 
more complicated. Indeed, as demonstrated in Fig. [12] 
where the fixed spin-down electron is successively placed 
away from the x-axis at (0, 20 nm) and at (0, -20 nm), the 
intrinsic structure corresponds to a superposition of two 
molecular isomers, each one described by a three-orbital 
wave function <&(§, |), but with the three localized spin- 
orbitals located on the vertices of two isosceles triangles, 
each one being a mirror reflection (relative to the x-axis) 
of the other. The base of the first isosceles triangle lies 
at -6 nm [Fig. 12 'a)] and that of the second one at 6 nm 



[Fig. 12 'a)] off the x-axis (in the y-direction) . 



The two-triangle configuration discussed for three elec- 
trons above may be seen as the embryonic precursor of 
a quasilinear structure of two intertwined "zig-zag" crys- 
talline chains. Such double zig-zag crystalline chains may 
also be related to the single zig-zag Wigner-crystal chains 
discussed recently in relation to spontaneous spin polar- 
ization in quantum wires j 32 l 33 l 

It is interesting to inquire of how this two-triangle 
structure is reflected in the spatial distribution of the 
electron densities. Indeed, in Fig. 13 'a), we display the 
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FIG. 14: (Color online) Spin-resolved conditional probability 
distributions for the (1/2,1/2) ground state of N = 3 electrons 
in an anisotropic dot at B — with parameters hu> x = 4.23 
meV, hujy = 5.84 meV (77 = 0.724), m* = 0.070m e and k = 1. 
(a) || CPD with the fixed spin-down electron located on the 
y-axis at (0,-20) (solid dot), (b) jl CPD with the fixed spin- 
down electron located off center at (40,11) (solid dot), (c) 
T| CPD with the fixed spin-down electron located on the x- 
axis at (43,0) (solid dot), (d) j j CPD with the fixed spin-ip 
electron located on the z-axis at (43,0) (solid dot). The spin of 
the fixed electron is denoted by a thick arrow (blue online). 
Lengths in nanometers. The vertical axes are in arbitrary 
units, but the scale is the same for all panels in this figure. 



electron density associated with the (3/2,1/2) ground 
state at B = 5 T. We note in particular the absence 
of a third peak at the center of the quantum dot. In- 
stead, two rather small peaks appear at (0,20 nm) and 
(0,-20 nm), in agreement with the two-triangle internal 
structure revealed by the CPD analysis. 



B. Strong repulsion (k = 1) 



We further display in Fig. 13 'b) the corresponding elec- 
tron density for the (1/2,1/2) ground state at B = 
and for a strong Coulomb repulsion (k = 1) at the in- 
termediate anisotropy r\ = 0.724. As a result of the 
enhanced electron localization, the electron density ex- 
hibits pronounced peaks whose locations form a clearly 
defined diamond; this indicates a gain the presence of a 
two-triangle internal configuration! 3 ^ The detailed inter- 
locking of the two triangular configurations is further 
revealed in the spin-resolved CPDs that are displayed 
in Fig. [TIJ From the CPDs in Figs, li^a) and [Tl^b), 



it can be concluded that one triangle is formed by the 
points Ri w (0,-20) nm, R 2 « (-43,10) nm, and 
R3 rs (43, 10) nm, while the second one (its mirror) is 
formed by the points Ri « (0, 20) nm, R' 2 rj (-43, -10) 
nm, and R^ « (43,-10) nm. The U [Fig. Qc)] and ft 
[Fig. [Hfd)] CPDS with the fixed electron on the right 
at (43^0) nm are similar to those in Figs.[8jb) and[8jc), 
respectively, with the difference that the central hump is 
clearly a double one. This indicates that each triangular 
configuration is associated with a wave function of the 
form $(A, |;1) given in Eq. (14 1. 

Naturally, the regime of a linear configuration versus 
a two-triangle one depends on both the strength of the 
interaction and the anisotropy. Detailed studies of the 
phase boundary between these two intrinsic structures 
are, however, left for a future investigation. 



VI. DEGREE OF ENTANGLEMENT 

The many-body wave functions for N — 3 electrons 
analyzed in the previous sections are highly entangled 
states, since they cannot be reduced to a single Slater 
determinant. For special ranges of the dot parameters, 
we showed that they acquire the same internal structure 
as the prototypical W-states. In this section, we demon- 
strate that the degree of entanglement can be further 
quantified through the use of the von Neumann entropy 
iS v n for indistinguis hable ferm ions which (in analogy to 
the two-electron cas e 35 * 36 * 37 * 38 ') is defined as 



S vN = -Tr(plog 2 p) + C, 



(17) 



where C is a constant (see below for choosing its value) 
and the single-particle density matrix is given by 



(V EXD \ala u \y EXD ) 

E^EXDjat^EXD)' 



(18) 



and is normalized to unity, i.e., Trp = 1. The Greek 
indices /i (or v) count the spin orbitals ip(v; p) that span 
the single-particle space (of dimension 2K; see Section 
|n| . Note that, in keeping with previous literature on two 
electrons] 35 * 37 * 38 ! the logarithms are taken to base two. 

Naturally, for calculating numerically the matrix ele- 
ments p vfll we use the expansion ^ to get 



EXD lata. 



I* 



EXD\ 



^^j(SD(/)|at a ,|SD(J)), 



i.j 



where the following conventions for the indices / (or J) 
apply: / - (ijk) and |SD(/)) = |^(1; i)V(2; j>(3; k)). 
The matrix elements (SD(/)|a| I aJSD( J)) between Slater 



determinants that enter in Eq. (19 1 simply equal ±1 or 
vanish. The single-particle density p in Eq. ( 18 1 is in 
general non-diagonal. Thus we further perform numeri- 
cally a diagonalization of p, and we use the new diagonal 
elements j5„„ to straightforwardly calculate the von Neu- 
mann entropy in Eq. (17). 



11 



As was discussed in Refs. 1361371381 the von Neumann 
entropy provides a natural measure of entanglement in 
the case of interacting indistinguishable fermions. In this 
case, the entanglement is related^ to quantum correla- 
tions that are intrinsic to the many-body wave function, 
i.e., £>vn quantifies the fact that strongly correlated states 
comprise a larger number of significant Slater determi- 
nants compared to weakly correlated ones. Accordingly, 
one expects that 5 v n increases when the many-body cor- 
relations increase (i.e., when R\y and B increase). This 
was the case indeed for the N — 2 quantum dot j 37 * 38 * 
but we have found that it also holds true for the N = 3 



1G 



quantum dot, as can be seen from Fig. 15 and Fig. 

In the case when the many-body wave function reduces 
to a single Slater determinant, i.e., when the expansion 
coefficients reduce to Ai = 5ij n , all the matrix elements 
vanish except three diagonal ones (corresponding to 
three fully occupied spin orbitals) which are equal to 
1/3; then — Tr(plog 2 p) = log 2 3 = 1.5850. Since the en- 
tanglement due to the Pauli exchange principle by itself 
cannot be used as a resource for quantum-in form ation 
processing] 36 ! 40 ! we take the constant C in Eq. (17) to be 



C= -log 2 N, 



(20) 



and as a result the von Neumann entropy for a single 
Slater determinant vanishes in our convention. 

In Fig. 1 1 5 1 we plot the von Neumann entropy for the 
three lowest EXD states with S z = 1/2 as a function of 
k (Rw) for N = 3 electrons in an anisotropic quantum 
dot for a strong anisotropy with anisotropy parameter 
X] = 1 /2. It is apparent that the von Neumann entropy in- 
creases for all three states as Rw increases (k decreases) 
and the electrons become more localized. 

At k = 12.50 (corresponding to weaker correlations), 
the von Neumann entropies for the three states are clearly 
non- vanishing, indicating that these EXD states are far 
from being close to a single Slater determinant. On 
the other hand, it is natural to expect that the EXD 
states will reduce to single Slater determinants at the 
non-interacting limit. To check this expectation, we have 
carried out an EXD calculation for the same QD param- 
eters described in the caption of Fig. |15[ but with a very 
large k — 10000 in order to approximately mimick the 
non-interacting limit. In this latter case, we found that 
indeed the ground state [with (S = 1/2; S z = 1/2) is 
practically a single Slater determinant made out from the 
three spin orbitals (to = 0,n = 0;|), (to = 0, n = 0; J.), 
and (to = l,n = 0; f) [the lowest-in-energy spatial or- 
bital (to = 0, n = 0) being doubly occupied; see Section 
[n] for the meaning of indices to and n). We also found 
that the corresponding 5 v n is practically zero. 

However, due to the r/ = 1/2 anisotropy, one has 
2Hlu x — hujy, which gives rise to a high degree of en- 
ergy degeneracy among excited Slater determinants with 
good total spin. For example, the Slater determinant 
| (to = 0,n = 0;T), (m = l,n = 0; J.), (to = l,n = 0; f)) is 
degenerate to the determinant |(m = 0, n = 0; |), (to = 
0, n = 0; I), (to = 0, n = 1; f)). In this situation, a small 
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FIG. 15: (Color online) Von Neumann entropy at zero mag- 
netic field for the three lowest EXD states with S z = 1/2 
as a function of the dielectric constant k [or equivalently 
the Wigner parameter Rw', see Eq. |S|] for N = 3 elec- 
trons in an anisotropic quantum dot with strong anisotropy 
(anisotropy parameter r\ = 1/2). Parameters: external con- 
finement huj x = 3.137 meV, hu> y — 6.274 meV; effective mass 
m* = 0.067ra e . The single labels 3/2 and 1/2 denote the 
quantum numbers for the total spin. The ground state (GS), 
and the first (I) and second (II) excited states are indicated. 
The horizontal arrow indicates the direction of increasing cor- 
relations. According to our convention, the von Neumann en- 
tropy for a single determinant vanishes. Although the energy 
gaps between the three EXD states diminish with decreasing 
k (they are quasidegenerate for k = 1), the relative energy 
ordering remains unchanged in the plotted range. 



e — e interaction is sufficient to produce strong mixing of 
the degenerate Slater determinants, and as a result the 
corresponding 5 v n values for excited states were found to 



be non- vanishing. These findings are reflected in Fig. 15 
where, for k = 12.50 (weakest Coulomb repulsion in the 
plotted range), the 5 v n value for the EXD ground state 
is noticeably lower than the values for the two excited 
states. 

As was demonstrated in Section |IV| at zero magnetic 
field and strong Coulomb repulsion (e.g., k = 1), the 
three electrons are well separated and localized, and their 
EXD wave functions are equivalent to the forms given in 
Eq. (03) (GS), Eq. ([XsJ (first excited state, I), and Eq. 
(13 1 (second excited state, II). These forms are special 
cases of the general form in Eq. ( 12 1 for which another 



measure of entanglement, called the tangle and specify- 
ing the reduced tripartite entanglement among the three 
localized spin qubitsp^ can be applied. 

The tangle can be calculated^ from the coefficients a, 
b, and c, and it was found that it vanishes for all cases 
covered by the general form in Eq. (12 1. In this respect, 



the von Neumann entropy for three well separated elec- 
trons studied here exhibits qualitatively a very different 
behavior, since the values of S* v n at k = 1 are all differ- 
ent, as seen from Fig. 15 In particular, we note that in 
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creases for all three states as the magnetic field increases 
and the electrons become more localized. An interesting 
feature is the discontinuous jump (around B — 3.4 T) 
in the von Neumann entropy of the EXD ground state. 
This jump is illustrated by the vertical arrows and is 
associated with the triple ABC point in Fig. [5] This dis- 
continuity in the ground-state SVn arises from the sudden 
change in the intrinsic structure (in term of constituent 
Slater determinants) of the ground state, as the latter 
changes its quantum numbers first from (1/2, 1/2; 1) to 
(1/2, 1/2; 2), and then again immediately to (3/2, 1/2) at 
the triple point. 

VII. SUMMARY 



FIG. 16: (Color online) Von Neumann entropy for the three 
lowest EXD states with S z = 1/2 as a function of the magnetic 
field for N = 3 electrons in an anisotropic quantum dot with 
strong anisotropy (anisotropy parameter r\ — 1/2). Parame- 
ters: external confinement hu> x — 3.137 meV, hu) y = 6.274 
meV; dielectric constant k = 12.5; effective mass m* = 
0.067m e . The single labels 3/2 and 1/2 denote the quan- 
tum numbers for the total spin. The vertical arrows indi- 
cate the discontinuous jump in the von Neumann entropy 
of the ground state at B = 3.4 T, where the ground-state 
quantum numbers change character, first from (1/2, 1/2; 1) to 
(1/2, 1/2; 2), and then immediately to (3/2, 1/2). The ground 
state (GS), and the first (I) and second (II) excited states are 
indicated, both to the left and to the right of the vertical ar- 
rows. According to our convention, the von Neumann entropy 
for a single determinant vanishes. 



this case the EXD value of SVn for the I state is lower 
than that of the GS state; this naturally reflects the fact 
that the first excited EXD state in this limit is effectively 
composed of only two Slater determinants [see Eq. (15)] 
compared to the three Slater determinants associated [see 



Eq. ( |14| )] with the EXD ground state. 

we plot the von Neumann entropy for the 



In Fig. 16 



three lowest EXD states with S z = 1/2 as a function of 
the magnetic field for N — 3 electrons in an anisotropic 
quantum dot with the same parameters as those for 
the energy spectra in Fig. [5] (strong anisotropy with 
anisotropy parameter r\= 1/2, as also was the case with 
Fig. 15 ). It is apparent that the von Neumann entropy in- 



We have presented extensive exact-diagonalization cal- 
culations for N = 3 electrons in anisotropic quantum 
dots, and for a broad range of anisotropics and strength 
of inter-electron repulsion. We have analyzed the exci- 
tation spectra both as a function of the magnetic field 
and as a function of increasing anisotropy. A main find- 
ing was the appearance of triple-crossing points in the 
ground-state energy curves for stronger anisotropics. 

Analysis of the intrinsic structure of the many-body 
wave functions through spin-resolved conditional proba- 
bility distributions revealed that for all examined cases 
(including those with parameters corresponding to cur- 
rently fabricated quantum dots) the electrons localize 
forming Wigncr molecules. For certain ranges of dot 
parameters (mainly at strong anisotropy), the Wigncr 
molecules acquire a linear geometry, and the associated 
wave functions with a spin projection S z — 1/2 are sim- 
ilar to the socalled IF-states that are a prototype of en- 
tangled states. For other ranges of parameters (mainly at 
intermediate anisotropy) , the Wigner molecules exhibit a 
more complex structure consisting of two mirror isosceles 
triangles. This latter structures can be considered as an 
embryonic unit of a zig-zag Wigner crystal in quantum 
wires. 

Finally, we demonstrated that the degree of entangle- 
ment in three-electron quantum dots can be quantified 
via the von Neumann entropy, in analogy with studies 
on two-electron quantum dots. 
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